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Context. The dependence of galaxy star formation activity on environment - especially in clusters - at high redshift is still poorly 
understood, as illustrated by the still limited number of z > 1.4 clusters on the one hand, and by the still debated star formation-density 
relation at high redshift on the other hand. 

Aims. The Zpimi ~ 2.2 JKCS 041 cluster allows to probe such environmental dependence of star formation activity at an unprecedented 
combination of redshifts and environments. Its study permits to enlarge the knowledge of high redshift clusters and to put strong 
leverage on observational constraints for galaxy evolution models. 

Methods. We analyze deep u'g'r'i'z'JHK, images from the CFHTLSAVIRDS surveys, which cover JKCS 041 cluster field. We 
first estimate photometric redshifts based on multi-wavelength photometry. We then lead a careful analysis to test the presence of a 
Butcher-Oemler effect. We work on galaxies within 2 x ^200 and with masses > 1.34 x 10"Mq, and use two comparison clusters at 
z = and z = 1 of similar mass. We estimate the radial profiles of the fraction of blue galaxies, taking into account the star aging with 
decreasing redshift. 

Results. After confirming the high redshift nature of JKCS 041, we find no evidence for a Butcher-Oemler elfect between z ~ 2.2 and 
z ~ for galaxies more massive than 1.34 x 10" Mq. In the cluster center, a change greater than Ayi;,„/Az = 0.16 between z ~ and 
z ~ 2.2 would be easily detected. We also find that JKCS 041 shows a consistent and systematic increase of the fraction of star-forming 
galaxies with cluster-centric distance, hence with decreasing density, for both a M > 1 .34 x 10" Mq selected sample and a lower mass 
sample. In particular, very few (less than 15%) star-forming galaxies are found within ?'2oo/2 among high mass (M > 1.34 x 10"Mo) 
galaxies. 

Conclusions. Our results show that the present-day star formation-density relation is already in place at z ~ 2.2. 
Key words. Galaxies: clusters: general - Galaxies: clusters: individual: JKCS 041 - Galaxies: evolution 



1. Introduction 

It is known that, in the local Universe, star formation activ- 
ity is correlated with galaxy environments: galaxies with low 
star formation rates are preferentially found in dense environ- 
ment, gakxiesjn cluster cores being virtually all quiescen t 
(e.g. IOemlel^[T974t iHogg et al.l I200I kauffmann et al.l |2004 . 
However the various processes which led to this s ituation in 
the lo cal Universe are still not fully understood (e.g. iTreu et all 
l2003h . To put constraints on how this star formation-density 
relation, well established in the local Universe, has been built 
through cosmic ages, one should observe galaxies at increasing 
look-back times and study their star formation activity with re- 
■ spect to the environment at fixed stellar mass, so as to isolate the 
' role of environment. In this context, galaxy clusters bring deci- 
sive informations, as they are the densest environments in the 
Universe. 

It has been shown t hat the star form ation-density relation 
holds out to z < 0.8 (e.g. lPatel et ani201lh . At z ~ 0.8-1 while 
studies in lo w-density environments (e.g. lElbaz et al 1 12007t 
ICooper et al.l i2008) observe a reversal of the star formation- 
dens ity relation, other studies focusing on cluster environments 
(e.g. iPatel et al.ll2009t iKovama et alll2010l) find that cluster core 
regions are devoid of star-forming galaxies. ISobral et al.l (1201 ll) . 
studying a very wide range of environments, nicely reconcile 
those observations: star formation activity increases with in- 
creasing density up to E ~ 10-30 Mpc"^, and then decreases 
with increasing density for 2 > 30 Mpc"^. 



Studying a superstructure at z ~ 1 .2, iTanaka et al.l (l2009h 
find that the star formation-density relation is already in place. 
When going to higher re dshifts (z ~ 1.4 - 1.6), t he sit uation is less 
clear On the one hand, iHavas hi et al.l (1201 Ol) and iHilton et"aD 
( 1201 Oi) study th e XMMXCS J2215.9-1738 cluster at z ^ 1.46 
{kT ~ 4.1 keV. IStanford""etani2006t [Hilton et al.ll2010l) : look- 
ing at the 24 //m and [Oil] emission respectively, both work ob- 
serve a high star formation activity at its centre. On the other 
hand, studies on the massi ve XMMU J2235-2557 cluster at 
z = 1.39 (/tr ~ 8.6 keV. |jee"et al.ll2009t iRosati et al.ll2009l) 
observ e that its core region does not present s t ar formation ac 
tivity 



(iLidman et al 
l20I0l:lBaueret al. " 



20081: iRosati et al l l2009t IStrazz ullo et jl 
:al.l(l20Ill) andlOuadri et jl 



- ■ 201 II) . IChuter et al.l (1201 1 

(120 1 II) study a large sample of galaxi es and environments i n the 
UKIDSS Ultra-Deep Survey (UDS, Lawrence et aDl2007l ; O. 
Almaini, in preparation) and find t hat the star fo r matio n-density 
relation holds out to z ~ 1.5-1.8. lOuadri et al] (|20I II) also in- 
vesti gate the CIG J() 2I8.3-05I0I cluster at z = 1.62 {kT ~ 1.7 
keV. lPapovich et al.ll20IOt ITanaka et al.ll2010l) and observe that 
its central region has an elevated fraction of quiesc ent objects 
relativ e to the field, in apparent disagreement with iTran et aH 
(120 1 Ol) . who observe an increasing of the star formation activ- 
ity along with the density. 

Though the two clusters XMMU J2235-2557 and XMMXCS 
J2215.9-1738 lie at a similar redshift of z ~ 1.4, XMMU J2235- 
2557 is massive and presents a well defined red sequence down 
to faint galaxies, thus being likely in a very advanced evolution- 
ary stage, whereas XMMXCS J2215. 9-1738 is less massive and 
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shows a deficit of faint galaxies on the red sequence, thus being 
in a less evolved dynamical state. At lower redshift [z ~ 0.3), 
[Braglia et al. (2009) studied the star formation activity in two 
clusters with opposite dynamical states: even if low star forma- 
tion activity is found in both cluster cores, a star formation activ- 
ity is found in the less evolved cluster, out to its virial radius and 
beyond, while no star formation activity is found in the more 
evolved cluster, thus hinting to a link between the dynamical 
state and the star formation activity in clusters. The different dy- 
namical status may explain the conflicting evidence observed at 
high redshift. 

The evolution of star formation acti vity within clusters 
has been studied in many papers (e.g. iHaines et al.l l2009l 
and references therein), starting with the pioneering work of 
iButcher & Oemled (11984 . The latter authors looked at the frac- 
tion of blue galaxies, fbiue, in clusters and its evolution with red- 
shift, thus probing the impact of dense environments on star for- 
mation activity. This seminal study found an increase of fbi„e 
in clusters with increasing redshifts - the so-called Butcher- 
Oemler effect - out to z ~ 0.5, thus pointing to an accelerated 
evolution in clusters. However, it has been subsequently shown 
that Butcher-Oemler eff'e ct studies may be sever ely aff'ected 
by methodological biases. ^Andreon & Ettoril ( 19991) have shown 
that a strong bias in the^utcher & O emleij ( 19841) cluster sample 
may account for the seen effbct. Acc ording to subsequent works 



dPe Propris et al.ll2004t Gotol 2005h. fbiue s eems not to depend 

is et al] Hi 



on the cluster mass. iDe Propris et al .1 (120031) have shown the ne- 
cessity to use mass-selected galaxy sample: in fact, if on e use an 
optica l luminosity-selected galaxy sample as Butcher & Oemled 
(Il984l) did, higher redshift samples will be biased towards low 
mass starburst galaxies, not included in lower re dshift samples, 
leading to an artificial increase of fbiue- At last, lAndreon et al.l 
(I2006h have shown that the criterium used to define a blue galaxy 
needs to take into account the younger mean age of the Universe 
and the secular increase in the star formation rate with redshift. 

Recent studies find no evidence for a Butcher-Oemler effect 
out to z ~ 0.5 ( Andreori et al, 2 006; Haines et al. 2009); nev- 



ertheless, lAndreon et al. (l2008bl) found evidence for a Butcher- 
Oemler effect when comparing a z ~ 1 cluster with local clus- 
ters. 

We here take advantage of CFHTLS/WIR DS public deep 
images covering the z ~ 2.2 JKCS041 cluster (lAndreon et alj 
|2009|) to address the afore mentioned issues. JKCS 041 presents a 
well defined red sequence (lAndreon & Huertas-Companvl201 ll) 
populated by a homogenous popu lation of galax ies with ex- 
tremely synchronized stellar ages ( lAndreonll201 1[) and an ex- 
tended X-ra y emission with T - 7.3^,g keV ( Andreon et alJ 
l2009ll20Tlh . attesting the presence of a formed potential well, 
deep enough to be hot and retain the intracluster medium. 
JKCS 041 thus offers a unique opportunity to probe star for- 
mation activity in clusters and the Butcher-Oemler effect out to 
z ~ 2.2. 

The plan of this paper is as follows: we describe in Section|2| 
the data used for JKCS 041, along with the analysis led on them. 
We estimate in Section [3] JKCS 041 photometric redshift. We 
then study the Butcher-Oemler effect in Section|4|and JKCS 041 
star formation activity in Section [5] We summarize and discuss 
our results in Section|6] 

In this paper, we adopt Ho - 70 km s ' Mpc ', Q„, = 0.30 
and Oa = 0.70. All magnitudes are in the AB system and masses 
are computed with a Chabrier (2003i) Initial Mass Function 
(IMF). JKCS 041 virial ra dius, estimated from the X-ray tem- 
perature, is r2oo = 1-53' dAndreon et alJl2009t ~ 0.76 Mpc at 



z ~ 2.2) and the cluster center is defined as the barycenter of the 
X-ray emission. 

2. Data and analysis 

2.1. Data 

JKCS 041 is in the ~ 0.6 deg^ area covered by CFHTLS deep 
survey (u'g'r'i'z' bands) and by WIRDS follow-up in the in- 
frared filters (JHKs bands, 50% point source completeness: 
Ks = 24.7) (Bielby et al., in preparation, catalogs are avail- 
able on the CFHT Science Data Archive siteQ). Throughout this 
work, we use the T0002 release of catalogs generated using 
band as detection image and the other bands in analysis mode. 
More specifically, we use magnitudes measured in 2" apertures 
(2" are ~ 17 kpc at z ~ 2.2) for colours, and "total" magni- 
tudes, both corrected for Galactic extinction using lSchlegel et akl 
(119981) . 

In order to use these catalogs for the study of JKCS 041, we 
need: a) to identify stars; b) to correct the underestimate of pho- 
tometric errors listed in the original catalog; c) to correct for 
(minor) residual photometric offset; d) to measure photometric 
redshifts and e) to correct for their systematic biases. We detail 
them in turn. 

The VIM OS VLT Deep Survey project (VVDS, 

iLe Fevre et al.l l2005h gives spectroscopic redshifts, Zspec, 
of several thousand of objects in the same area. For our spectro- 
scopic sample, we use objects in common with T0002 catalogs 
- rejecting the edgings - with < 23 (see hereafter) and a 
secured Zspec (flag=3,4), thus yielding a spectroscopic sample of 
2537 galaxies {Zspec ^ 1-5) and 366 stars {Zspec = 0). 

2.2. Star removal 

St ars are identified, an d removed, in the colour-colour plane, as 
in ICowie et al ] (11994 and later works. In the z' - vs g' - z' 
colour-colour diagram, spectroscopic-identified stars populate a 
narrow locus offset from galaxies, as shown in Figure [Tj We 
therefore classify as star every object bluer in z' - than the 
broken line shown in Figure [1] This criterium excludes more 
than 95% of the spectroscopically confirmed stars. We checked 
that the latter constitute an unbiased (for our purposes) sample 
of stars, because they cover the sa me colour-colo ur locii of the 
large and representative sample in iFinlator et al.l (HoOfl). Stars 
not identified as such in this phase will be removed in later 
phases (during the photometric redshift selection and the back- 
ground statistical subtraction phase). 



2.3. Photometric errors correction 

Because of a slight noise correlation introduced by image 
resampling during the stack i ng, fl ux errors as derived by 



SExtractor dBertin & ArnoutsI 1996|) ar e underestim ated (e.g. 
ICa sertano et all 120001: lAndreon! 1200 ll) . Following lAndreonI 
(200jr by binning the images we can recover the actual back- 
ground noise, and thus compute the underestimation factor, 
thanks to the fact that the correlation is present on small spa- 
tial scales only . We found a facto r o f 1.5 for optical band s, in 
agreement with lllbert et al.l(l2006l) and lCoupon et all (l2009l) . and 
a factor 2.0 for near-infrared bands. The larger factor for near- 
infrared bands is due to the native larger pixel size of the images 
and thus more finely resampHng in the stacking. After correction 



,http : //wwwl ■ cade- ccda . hia- iha . nrc- cnrc . gc . ca/ cfht/WIRDSTQC 
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Fig. 1. Star/galaxy separation: g' - z' vs z' - Ks diagram. 
Noisy contours give the density of objects at each location. 
Spectroscopically confirmed stars are marked with black dots. 
The broken line marks the star/galaxy separation threshold. 

Table 1. Systematic offsets nifi, - m,neas between measured and 
best-fit model magnitudes 



S 



J 



H 



-0.08 0.06 0.01 0.03 -0.01 -0.05 -0.02 



-0.01 



of the error underestimate, completeness, defined as S /N - 5, 
occurs at Kg - 23 mag, J = 23.6 mag andz'-7 = 29.5-1.25x/r.5 
mag. Lower S /N data are never used in this work and even more 
restrictive cuts are used in most instances, as detailed below. 



2.4. Photometric redshift estimation 

For photometric redshift estimation, we use Eazy 
dBra mmer et al. 2008) with default settings and a Ks-hand 
magnitude prior. We use both the full photometric redshift 
probability distribution function, p{z), and Zmp, the latter as 
a point estimate of the photometri c redshift. Zmo is the red 



shift posterior mean (see §2. 5 of iBrammer et al.l l2008l). As 
noted in previous works (e.g., Brodwin et al I 2006t Ilbert et ap 
l2006l ICouporretaDl2009l lllbert et al.ll2009t iBarro et al.ll201 Ih . 



off'sets in the photometric calibration or the inadequacy of 
the templates to reproduce all the observed SEDs can lead 
to systematic offsets in the photometric redshift estimation. 
Eazy addresses the latter by using a template error function. 
In order to fix photometric offsets, we use the subsample of 
our WDS galaxies having a S /N greater than 10 in all the 8 
photometric bands (2281 galaxies) and we compute the mean 
difference between the predicted (best fit model) and observed 
magnitude. More precisely, we compute this average one filter 
at a time, we apply the photometric offsets, and iterate until 
the procedure converged. Five iterations were sufficient for 
convergence and the found offsets are listed in Table [T| These 
small shifts, comparable with previous works (lllbert et al.ll2006i 
l2009t ICoupon et all 120091: iBarro et al.ll20Tll) . were applied for 
photometric redshift estimation only. 

Figure |2] shows Zmp vs Zspec (top panel) after coiTection of 
the photometric offsets. The scatter is 0.019 + 0.098, much bet- 
ter than if the correction is not applied (0.030 + 0.097), largely 
because an improvement at low redshift, as shown from com- 
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Fig. 2. Performance of the photometric redshift estimate: con- 
tours give the density of objects at each location. In this Figure 
only, we consider that the photometric redshift estimate fails 
when Zspec does not belong to the 3cr confidence interval of Zmp 
(light magenta). Upper panel: Zspec vs Zmp when offsets are ap- 
plied. Middle panel: residuals A = {Zspec - Zmp)/{^ + Zmp) vs 
Zmp when offsets are applied. The light blue line is a linear fit 
to the three points in light blue, representing the median values 
when the data are binned in three bins. We observe that Zspec S 1 
galaxies tend to have their Zmp value underestimated (cf. light 
blue line). Lower panel: residuals A vs Zmp when no offsets are 
appHed. 



parison of residuals before and after coiTection (see middle and 
bottom panels). 

Even after correction of the photometric offsets, the photo- 
metric redshift tends to underestimate the spectr oscopic redshift 
at z > 1 (see middle panel), as already noted by iBrammer et al.l 
(l2008l) . The linear fit to the binned data (light blue) illustrates 
this trend. We checked that a similar underestimate of the red- 
shift holds for the sub-sample of red galaxies, defined as {U - 



vy 



model 
rest— frame 



> 1.3 mag. We apply this last correction in the only 



place of this work where it is needed, in Section |3] 



2.5. Background removal 

When estimating the properties of JKCS 041 galaxy population 
- in Sections |4| and jSj we need to account for galaxies on the 
cluster line of sight, that we generally call background, by a two 



3 



A. Raichoor & S. Andreon: Star formatio 

Step procedure: first, we will make a photometric redshift se- 
lection by removing galaxies (and QSO) that are at z < 1.7 or 
z > 3.5 at > 99 % confidence. This selection is accomplished by 
keeping objects with: 

Jr»1.7 r»+oo 
p{z)dz< 0.99 and ( piz)dz< 0.99 (1) 
J3.5 

We note that the latter equation removes only very few objects. 
We remark that this selection is very effective in removing "low 
redshift" galaxies, as it removes 2495 of the 2522 galaxies with 
Zspec ^ 1.5 in our VVDS sample. As later detailed, we explored 
other possible choices, and results are insensitive to the precise 
used recipe. 

As a second step of background subtraction, we use a large 
control area (~ 0.1 deg^ around the cluster, excluding a disk of 
7' radius centered on the cluster) to estimate the residual back- 
ground. This step also subtracts any star not identified as such 
by colours. 

We here analyze the impact on our sample selection of wrong 
photometric redshift estimate. In Figure |2] a sizable number of 
points scatters off from the diagonal in the top panel. In this 
Figure only, we consider that the photometric redshift estimate 
fails when Zspec does not belong to the 3(t confidence interval 
of Zmp (galaxies in light magenta). A large majority of the out- 
liers are objects with fairly large errors. Most of the outliers 
are in the lower-right corner of the Zspec vs Zmp plot (i.e. galax- 
ies with overestimated photometric redshift). These galaxies in- 
crease the noisiness of our measurements but do not introduce 
any bias, because kept in the sample. The most troublesome ones 
are high-redshift galaxies with a largely underestimated redshift 
(i.e. the objects in the very top-left corner) and a nominal small 
redshift uncertainty. In fact, these galaxies, if they exist, would 
be a source of incompleteness in our sample, which discards all 
galaxies which are at z < 1.7 at 99% confidence. In our spectro- 
scopic sample, this very situation never happens; there are only 
two galaxies - lying at a redshift {Zspec ~ 1 -6) lower than the one 
we are interested in - which are in a qualitative similar situation. 

3. JKCS041 photometric redshift estimate 

In order to estimate JKCS041 redshift, we select bright {K, < 
21.2) objects within 0.5 x r2oo and within 3cr from the z - J 
vs K, colour-magnitude relation (1.74 < z - J < 2.2. lAndreonl 
I2OI lb . because this choice maximizes the cluster membership 
likelihood (such bright and red galaxies are rare in the field, as 
measured all around the cluster). We emphasize that we do not 
use any photometric redshift selection in this Section. For each 
of these 8 galaxies. Figure [3] shows their SED along with the 
Eazy best-fit template and, as insets, its position in the colour- 
magnitude diagram (upper-left inset) and the photometric red- 
shift probability distribution function p(z) [lower-right inset). 
We observe that the fits are of good quality and that those 8 
galaxies show a prominent 4000 A break near the J band, char- 
acteristic of high-redshift old galaxies. We remark that the two 
bluest galaxies have a p(z) less peaked and slightly shifted to- 
wards lower redshifts. 

In Figure H] we gather the photometric redshift probabil- 
ity distribution functions p{z) for the 8 selected galaxies (thin 
coloured lines). All the 8 p{z) peak between z ~ 1.8 and 
z ~ 2.2. If we assume that all those 8 objects belong to the 
cluster, the cluster photometric redshift is obtained by multi- 
plying the p{z) functions, as they have been derived from inde- 
pendent data. The result is plotted in Figure |4] as a thick black 
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Fig. 4. JKCS 041 photometric redshift estimate : the thick black 
line shows the cluster photometric redshift, the thin coloured 
lines represent photometric redshifts of the 8 bright galaxies on 
the cluster red sequence shown in Figure [3] The area shaded in 
yellow represents the Icr confidence interval. The top axis shows 
the raw photometric redshift, not accounting for the known un- 
derestimate of photometric redshifts. 



line: it has a very peaked shape around z ~ 2.0 (Icr confi- 
dence interval: [1.97,2.02]) in broad agreement with the pho- 
tometric redshift of z„hot - 2 - 20 + .11, robustly estimated in 
lAndreon & Huertas-Companvl (1201 ih . We stress that our pho- 
tometric redshift estimate is derived using a photometric red- 
shift correction extrapolated from z < 1.5: it i s henc e slightly 
less robust than lAndreon & Huertas-Companvl (1201 ll) estimate, 
which uses the red sequence and is calibrated on the cluster CIG 
J0218.3-0510, which has a spectroscopic redshift Zspec - 1 -62 
(iPapovich et al.ll2O10HTkiaka et al.ll2010l) . 

The top axis of Figure |4] shows that JKCS 041 would have a 
high redshift {Zrmv ~ 1 -8) even when neglecting the correction 
for the underestimate of photometric redshift shown in Figure |2] 
This underestimate hints to a calibration problem for the models. 
Nevertheless, as model predictions match JKCS 041 colours for 
a redshift z ~ 1 .8, we will use in the rest of this article this value 
to redshift our models. This approach allows to use predicted 
colours in agreement with our observations (in particular, see 
Figure|9]where predicted colours for a Simple Stellar Population 
match the red sequence). 



4. Butcher-Oemler effect 

The measurement of the Butcher-Oemler effect requires partic- 
ular attention to the way galaxies are selected in mass and clas- 
sified as blue/red (see Section [T]i. First, we use a mass-selected 
sample. We select all galaxies more massive than 1.34x 1O"M0, 
which corresponds to our completeness {S /N = 5) limit in the 
worst case (old red galaxies). By mass we refer, as in previous 
works, to the Bruzual & Chariot (2003)) models (2007 version, 
CB07 hereafter) model mass, and specifica lly the mass of the gas 
that w ill eventually turn into stars (see e.g. lLonghetti & Saraccd 
(I2OO9I) for other possible definitions of model mass). These 
(model) masses are computed for solar metallicity, formation 
redshift z/ = 5, and either Simple Stellar Populations or expo- 
nentially declining star-forming t model with < SFH t (Gyr) 
< 10. 
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Fig. 3. JKCS041 bright red sequence galaxies analysis. Large panels shows the normalized SED (black dots with error bars), 
along with the Eazy best-fit template (cyan line) and the best-fit model photometric points (red dots). The thin light gray shaded 
area around the best-fit template shows the model uncertainty in the best-fit template, as given in Eazy. Upper-left insets: Colour- 
magnitude relation of all galaxies within 0.5 x rjoo- The closed (red) point emphasizes the object whose SED is shown. Red dashed 
lines represent cuts in magnitude and colour used in this plot, whereas the black vertical dotted line represents the - 23 limit. 
Lower-right insets: photometric redshift probability distribution function p{z) for the considered galaxy. All those red and bright 
galaxies have a SED with a prominent 4000 A break and are well fitted with a template at z ~ 2.0. 



4. 1 . Blue/red definition 

We define a galaxy as blue if it is b l uer than a CB07 m odel 
with T ^ 3.7 Gy r, as in lAndreon et all ( |2004 l2006l l2008bh and 
iLoh et all (l2008h . This galaxy will be bluer by 0.2 mag mB-V 
than red sequence galaxies at the present e poch (which would 
be a b lue galaxy by the original definition of iButcher & OemleJ 
1984. The rationale behind this choice is described in detail in 
Andreon et al. I (l2006l) and lAndreonI (l2006l ). 3ut in short we take 
into account the stellar evolution of galaxies with time. 

We emphasize that a galaxy to be classified as blue should be 
very blue, at least 0.8 mag bluer mJ-Ks than the red sequence. 



This is imposed by the requirement of following galaxies back 
in time, with a criterium independent of the redshift (under the 
assumption of exponentially declining SFH). A significantly nar- 
rower selection is precluded because, at lower redshift, it would 
end up on the top of the red sequence. 

Figure|5]shows the / - A'j vs Ks (about rest-frame u* - r' vs. 
r') colour-magnitude diagram for galaxies within 2 x r2oo of the 
cluster center and more massive than 1.34x IO^Mq. We adopted 
this colour index because it best matches the colour index used 
in our comparison sample below. As later detailed, our results do 
not depend on the used colour index. Not many galaxies are blue 
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Table 2. Cluster sample used to estimate the Butcher-Oemler 

effect 



Fig. 5. Colour-magnitude diagram for galaxies in JKCS 041 line 
of sight (i.e. JKCS 041 anc/ residual background galaxies) within 
2 X r2oo- The black dotted lines represent the S /N - 5 com- 
pleteness. The green thick solid line represents the mass cut 
(M - 1.34 X 10"Mq). The dashed line indicates the blue/red 
threshold. 



(below the dashed line in Figure |5]l, even before accounting for 
residual (after photometric redshift selection) background in the 
JKCS 041 line of sight. In addition, we underline that all (resp. 
more than 90% of) our sample galaxies have a S /N > 10 in /T., 
(resp. J - Ks). 

4.2. Measurement of a Butcher-Oemler effect 

We compute the blue fraction and the number of red and blue 
members galaxies accounting for residual background galaxies 
(i.e. in the line of sight, and not belonging to the cluster) us- 
in g our control fie l d follo wing the Bayesian methods introduced 
in I Andreon et al.l (l2006h . We adopt uniform priors for the pa- 
rameters. We consider three regions, defined by r/r2oo ^ 0.5 , 
0.5 < r/r2oo < 1 and 1 < r/raoo < 2. 

In order to provide lower redshift reference clusters, w e 
consider RzCS052 at z = 1. 016 (lAndreon et al.l HoOSbh . 
and Abell 496 (A496 hereafter, lAbelll Il958h at z ^ 0.033 
dStruble & RoodI Il999l). already used in the Butcher-Oemler 
studv in lAndreon et al.l (l2008bh . considering a lower mass 
threshold. Table |2] lists some key characteristics of our cluster 
sample used to study the Butcher-Oemler eff'ect. As one can 
see, our three clusters have roughly similar masses; moreover 
our study shows that they also have similar richness (see here- 
after. Figure |6|l. We use as colour index u* - r for A496 and 
I - z' for RzCS 052, thus adopting in both cases a rest-frame 
M*-like band as blue band, as for JKCS 041. We recompute the 
blue fraction of A496 and RzCS 052 clusters for galaxies with 
M > 1.34 X 10"Mo. 

Figure |6] shows red (left panels) and blue (right panels) ra- 
dial number profiles (upper panels) and radial density profiles 
(lower panels) for all three clusters, after accounting for the 
background. The three clusters have similar blue and red radial 
profiles, showing that the three clusters have similar richnesses. 
The right panel shows that all these three clusters have a negli- 
gible number of blue galaxies within 2 x r2oo- Unfortunately, the 
low redshift of A496 prevents us to probe its outer regions. 

Figure|7]shows the (posterior) probability distribution of the 
blue fraction of the three clusters in the three radial regions (two 
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Fig. 6. Radial number profile (upper panels) and radial number 
density profile (lower panels) of red (left panels) and blue (right 
panels) member galaxies as a function of the cluster-centric dis- 
tance, coded as indicated in the legend, after accounting for the 
background. The three radial ranges used are indicated by the 
magenta lines. Points indicate the maximum a posteriori, the er- 
ror bars the shortest 68% interval. Points are sUghtly shifted hor- 
izontally for readability. 



for A496) and list their summaries (point estimates and 68 % 
shortest interval). We emphasize that our fbiue computation re- 
quires the use of the full probability distributions for intervening 
quantities, «i/„(. and iired, not just their point estimates and 68% 
uncertainties reported in Figure |6] 

Figure [8] shows the blue fraction profiles of the three clusters 
and our main result of this section: we observe the same (negli- 
gible) amount of blue galaxies more massive than 1.34x 1O"M0 
all the way up to z ~ 2.2 in all radial bins. The values of fbuw are 
all less than 1 .Act away for all the three clusters in all the radial 
bins. 

To quantify the sensitivity of our data to an evolving fbute, 
we focus on the innermost bin {r/r2oo < 0.5) of the two clus- 
ters at the redshift ends, A496 and JKCS 041, to maximize 
the redshift leverage. By looking at the posterior of AfhUie = 
fhiff""'' - fbiT^ ^fi'i'<e < 0.36 (resp. A/w„, < 0.18) with 95% 
(resp. 68%) probability. Hence the rate at which the blue frac- 
tion changes, Afhi„e/Az, is less than 0.16 (resp. 0.08) per unit 
redshift at 95% (resp. 68%) probability. By a way of compari- 
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Fig. 7. Posterior probability distribution for the blue fraction for 
the three clusters (JKCS 041: upper row, RzCS 052: middle row, 
A496: lower row), in three radial ranges. The 68% shortest con- 
fidence intervals are shaded in yellow. 
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Fig. 8. Blue fraction as a function of the cluster-centric distance: 
JKCS 041 (z ~ 2, black dots), RzCS 052 (z ~ 1, green stars), 
and A496 (z ~ 0, red squares). Points indicate the maximum 
a posteriori, the error bars the shortest 68% interval. The three 
radial ranges used are indicated by the magenta lines. Abscissae 
are slightly shifted for readability. 



son, iButcher & Oemled(ll984 found a slope of 0.5 for a roughly 
similar radial aperture, but for a sample including less massive 
galaxies. Although error bars at z ~ 2.2 are large, the redshift 
leverage of this work allows us to tightly constrain the evolution 

of f blue- 

4.3. Robustness of the result on assumption 

We test the robustness of our result on assumptions: a) we 
change the redshift of JKCS 041 by +0.1; b) in place of the 
selection given in Equation [T] we keep in the saniple g alaxies 
with at least 30 % probability (following fTran et al.ll20Ta) in the 
1.5 < z < 2.1 range, i.e. 



-2.1 



r piz)dz > 0.3 

Jl.5 



(2) 



Fig. 9. Colour-magnitude diagram for galaxies in JKCS 041 line 
of sight (i.e. JKCS 041 and residual background galaxies) within 
2 X r2(K). The black dotted lines represent the S /N = 5 locii. 
The green solid line represents the locus of constant mass M = 
1.34 X 1O"M0. The black dashed line indicates the locus of a 
constant SFR of 4 Mq yr"', splitting quiescent and star-forming 
galaxies. 



c) we use the z' - J (rest-frame -2500-3000 A - u") instead 
of J - Ks colour index to identify red and blue galaxies. For 
those three cases, none of our measurements (fractions and radial 
profiles) change, not even by Icr. If instead of a mass-selected 
sample, we use a luminosity-selected sample brighter than an 
evolved My - -20.8 mag for all three clusters, then we derive 
consistent radial profile for the three clusters, i.e. we continue to 
see no Butcher-Oemler effect. 



5. Star-formation activity 

In this section we classify galaxies as star-forming, or quies- 
cent, from the slope of the UV continuum, specifically according 
whether they are bluer or redder in thez'-7(~ ^2800/^^3700) vs A'^ 
diagram than a CB07 model of solar metallicity, z /om = 5 and 
an exponential declining star formation history with star forma- 
tion rate (SFR) equals to 4 Mq yr ' at z - 1.8. Th i s SFR value 
is chosen according to the analysis of iKriek et al.l (l2009h . who 
studied a z ~ 2 spectrum of a quiescent galaxy and found a max- 
imum SFR of 4 Mq yr ' (for a Chabrier IMF). We underline that 
there is no reference SFR value at z ~ 2 used to classif y a galaxy 
as quiescent: we chose the value of lKriek et al.l (l2009l) study, be- 
cause it relies on a spectro scopic measuremen t, though only on 
one object. Other works, as lOuadri et al.l (1201 ll) for instance, use 
a different criterium. In addition to our M - 1 .34x 1O"M0 mass- 
selected sample, we consider less massive galaxies inaS /N se- 
lected sample. The latter choice has the advantage of enlarging 
the mass range, but the disadvantage of making difficult, not to 
say impossible, to compare results derived from data of diff'er- 
ent depths or for clusters at diff'erent redshift. Figure|9]shows the 
data and the various relevant locii. 

Using the same Bayesian methods used in previous section 
for red and blue galaxies, we compute the radial profile of qui- 
escent and star forming galaxies, as well as the radial profile 
of the star-forming fraction, shown in Figure \W\ This Figure 
shows a consistent and systematic increase of the fraction of 
star-forming galaxies with cluster-centric distance for both the 
M > 1.34 X 10"Mo sample (left panels) and the sample of less 
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Fig. 10. JKCS 041 radial profiles of cluster quiescent rig {upper 
row) and star-forming «,/ (middle row) galaxies, and the cluster 
star-forming fraction /,/ (lower row) for the high mass subsam- 
ple (left column) and the less massive subsample (right column). 
The three radial ranges used are indicated by the magenta lines. 



massive galaxies (right panels). This is our main result of this 
section. 

Figure [TT| shows that local density, estimated from the dis- 
tance of the 7''' nearest neighbour, goes as a function of cluster- 
centric distance: it decreases with r/r2oo until r/r2()o ~ 1, af- 
ter which its trend cannot be estimated because of uncertainties. 
Note that the cluster center is defined by the X-ray barycenter, 
not by the peak of the local density itself and thus the peak at low 
radii is not due to a selection effect. Since density and cluster- 
centric distance run hand in hand, at z ~ 2.2, the fraction of 
star-forming galaxies decreases with density. For illustrative pur- 
pose, we plot in Figure[T2]JKCS 041 star-forming fraction fsf vs 
background-subtracted density Z7. 

We emphasize that while the high mass sample is complete 
(in mass), the sample of lower mass galaxies has a completeness 
depending on whether the galaxy is quiescent or star-forming. 
This prevents any quantitative comparison of the star-forming 
fraction values between the two samples, because, for exam- 
ple, a larger fraction of star-forming galaxies among less mas- 
sive galaxies may be genuine or just a selection effect, due to 
mass-incompleteness of quiescent galaxies of low mass, as al- 
ready mentioned. On the other hand, this selection effect is inde- 
pendent on the cluster-centric distance and thus does not affect 
our conclusion about the increase of the fraction of star-forming 
galaxies with cluster-centric distance. 

We test the robustness of our result on assumptions: a) we 
change the redshift of JKCS 041 by +0.1; b) we use Eql2]for 
photometric redshift pre-selection; c) we multiply/divide by 2 
the SFR threshold value used to define quiescent/star-forming 
galaxies; d) we classify galaxies in star-forming, or quiescent, 
in the U - V vs V - J plane, as iWilliams_et al . (200S) and 
lOuadri et al.l ( 1201 ll) . We emphasize that in case d) the two classi- 
fications are almost identical, but that our rest-frame J photome- 
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Fig. 11. Local density, estimated from the distance of the 7''' 
nearest neighbor vs JKCS 041 cluster-centric distance: for this 
Figure only, galaxies are selected solely with a photometric red- 
shift criterium (Eq. [1) and K s < 23 mag. The re d line shows 
the best-fit with a theoretical iNavarro et alJ (Il996h profile plus 
a constant to take into account the background. The three radial 
ranges used for estimating ftiue and /j/ are indicated by the ma- 
genta fines. 
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Fig. 12. JKCS 041 star-forming fraction fsf vs density for the 
high mass subsample (left column) and the less massive subsam- 
ple (right column). The plotted values for £7 are the mean values 
of the best-fit profile in Figure[TT](red line) for each radial range, 
after background subtraction. 

try comes from extrapolating the available multicolour photom- 
etry. Our result does not change with case a). For cases b), c) 
and d), we observe a radial profile increase with cluster-centric 
distance, i.e. that the star formation-density relation is akeady in 
place at z ~ 2.2. We stress that in all cases, the fraction of mas- 
sive star-forming galaxies within r2oo/2 is very low (< 15%, Icr 
errorbars included). 

6. Discussion and conclusions 

In this work, we took advantage of CFHTLS (u*g'r'i'z' bands) 
and WIRDS (JHR^ bands) images to study the JKCS 041 clus- 
ter. Using two lower redshift clusters of similar mass (RzCS 052 
at z - 1 .016 and A496 at z = 0.032) as a comparison sample, we 
studied the evolution with redshift of cluster galaxies properties. 
Our results are: 

1. By using photometric redshifts, we confirm that 
JKCS 041 is a high redshift cluster, in agreement with 
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lAndreon & Huertas-Companvl (1201 ih . Our estimate is 
Zphot = 2.00;^Q Q3 (Icr errors), after a systematic correction of 
6z = 0.2 extrapolated from z < 1.5. 

2. Working on a mass-limited sample (M > 1.34x 10"Mq) and 
taking into account the star aging with decreasing redshift, 
we observe the same (negligible) fraction of blue galaxies 
all the way up to z ~ 2.2 in all radial bins within 2 x r/r2oo ■ 
we thus do not observe any evidence of any Butcher-Oemler 
effect between z ~ 2.2 and z ~ 0. With the forced choice of 
our definition of blue galaxy, very few galaxies more mas- 
sive than M - 1.34 x IO^Mq are blue, for all redshifts and 
radii. Although error bars are large, the redshift leverage of 
this work is at least twice larger than any previous work, al- 
lowing us to reject with confidence a change greater than 
A/wac/Az = 0.16 in the cluster center. 

3. JKCS 041 shows a consistent and systematic increase of 
the fraction of star-forming galaxies with cluster-centric dis- 
tance for both the M > 1.34 x 10" M© sample and a sample 
of less massive galaxies. In particular, very few (less than 
15%) star-forming galaxies are found within r2oo/2 among 
high mass galaxies. 

4. Density decreases with increasing cluster-centric distance 
for r < r2oo, and thus our statements above may be rephrased 
in terms of local density. 

lAndreon et all ( l2008bl) led a very similar analysis of the 
Butcher-Oemler effect for RzCS 052 and A496, the only differ- 
ence being that a lower mass cut is used (M > 4 x IO'^Mq). 
Interestingly, they find an evolution in the fraction of blue 
galaxies between z ~ 1 and z ~ 0, whereas we do not. 
Gathering the ir and our results h i nts towards a downsizin^-like 
scenario (e.g. Cowie et al.lFl996l: iTreu et all 120051 : llovino et afl 
l2010HPeng et al.ll2010h . where the most massive galaxies have 
their properties set in the very early Universe (z » 1), while less 
massive galaxies still evolve in < z < 1 . 

Our observations of JKCS 041 show that most of bright core 
cluster galaxies are red and passive at z ~ 2.2, an observation that 
models still have difficulties to match (e.g. iMenci et al J [20081 : 
iRomeo et al ] l2008h . 

The presence of a star formation- density relat i on in 
JK CS 041 is in agreemen t with the works of lChuter etall (1201 lb 
and lOuadri et alJ (l201 ih at z < 1.8 and those on the XMMU 
J2235-2557 cluster at z = 1. 39 (Lidman e t al. 2008: Rosati et al.i 
120091: IStrazzullo et alj|20Tol: iBauer etldlbOl II) and in disagree- 
ment wit h studies on the XMMXCS J2215. 9- 1738 cluster at 
z= 1 .46 (iHavashi et al.ll2010t [Hilton et al.ll2010l) . 

This variety of results may be either a manifestation of 
a spread of star formation-density relations at high redshift - 
perhaps due to the cluster dynamical status, or just the result 
of unidentified systematics. Enlarging the sample used for the 
Butcher-Oemler effect - and particularly selecting likely low 
redshift descendants of JKCS 041 - and deepening the analy- 
sis of the star formation activity in JKCS 041 is necessary to 
strengthen our results, and will be addressed in future works. 
However, JKCS041 is uniquely adapted for galaxy evolution- 
ary studies, even if it has a photometric redshift only, because 
it is so far the only z > 1 .5 cluster with measured intra-cluster 
medium properties - hence rooo - and a sizable red population. 
More clusters, with robust estimate of mass are needed to con- 
solidate the link, put forth by JKCS041, between star formation 
and environment at high redshift. 
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